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Abstract 



We study a resistively shunted semiconductor superlattice subject to a high- 
frequency electric field. Using a balance equation approach that incorporates 
the influence of the electric circuit, we determine numerically a range of am- 
plitude and frequency of the ac field for which a dc bias and current are gen- 
erated spontaneously and show that this region is likely accessible to current 
experiments. Our simulations reveal that the Bloch frequency corresponding 
to the spontaneous dc bias is approximately an integer multiple of the ac field 
frequency. 
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In 1928 Bloch demonstrated that electrons in a periodic lattice potential with period 
a subject to a dc electric field E dc undergo oscillations with characteristic frequency ujb = 
eaEd c /h. In naturally occurring solids, observations of Bloch oscillations are precluded by 
the requirements of the extremely high applied electric fields needed to reduce the Bloch 
oscillation period below the dephasing times arising from the ever-present scattering due to 
impurities and phonons. In 1970, however, Esaki and Tsu || recognized that the inherently 
longer spatial periods available in semiconductor superlattices (SSLs) allow Bloch oscillations 
in these artificial structures to lie the the terahertz (THz) domain — comparable to or higher 
than the corresponding scattering frequencies — even for modest field strengths (~1 kV/cm). 
Thus, not only should the Bloch oscillations be observable in SSLs, but SSLs could also serve 
as devices to produce THz radiation 0. 

Since 1970, an enormous literature has evolved describing novel physical phenomena in 
SSLs and their potential applications for devices ||. Among the most recent developments, 
which have relied on advances in both electromagnetic radiation sources and coupling tech- 
niques, are detailed observations of the influence of a THz field on the dc conductivity of an 
SSL including (1) THz multiphoton-assisted tunneling [|[; (2) ac-field-induced reduc- 

tion of the dc current ||; (3) absolute negative conductance ||; and (4) resonant changes in 
conductivity [0]. 

In the present article, we examine a question that is roughly the converse of the Bloch 
oscillations in an SSL: namely, can a purely ac external field applied to an appropriately 
configured (but unbiased) SSL create a dc bias and corresponding dc current? In technical 
terms, this corresponds to a spontaneous breaking of spatial reflection symmetry. We show 
that the answer to this question is affirmative. Using the methods of nonlinear dynamics, 
we find the condition for spontaneous dc current generation in the circuit consisting of the 
SSL shunted by an external resistance. We observe that the Bloch frequency corresponding 
to the spontaneously generated bias across the SSL is approximately an integer multiple of 
the frequency of the external ac field. Finally, we use the time-periodicity of the applied 
field to understand the origin of this "phase- locking" . 

We consider an SSL of spatial period a, length I and area S that is homogeneously doped 
with N carriers per unit volume. The SSL is shunted by an external measuring device of 
resistance R and is exposed to an external ac electric field E ext (t) = E cos fit. We use the 
standard tight-binding dispersion relation for electrons belonging to a single miniband, so 
that e(p) = A/2[l — cos(pa/h)], with e being the electron energy, p the quasimomentum 
along the SSL's axis, and A the miniband width. We describe the dynamics of electrons by 
the balance equations |8]|| 

V = -eE tot (t)/m(S) - 7t ,V, 8 = -eE tot (t)V - le (8 - £<®), E sc = 4-7reNV/e - aE sc 

(1) 

where V(t) and 8(t) are the electron velocity and energy averaged over the time-dependent 
distribution function satisfying the Boltzmann equation, and 7^ and 7 e are phenomenological 
relaxation constants for the average velocity and energy, respectively. Elastic scattering by 
impurities, interface roughness and structural disorder contributes to j v ; inelastic phonon 
scattering comprises the main channel of energy dissipation. The parameter 8^ describes 
the electron equilibrium energy ||, which is a function of the lattice temperature. Also 
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appearing in (0) is the energy- dependent effective mass m(8) = mo(l — 28/ A) x , where 

<m — 27i2 

The electric field acting on the electrons, E tot {t) = E sc {t) + E ext (t), is the sum of the 
external ac field E ext (t) and the self-consistent field E sc (t), which incorporates the influence of 
the circuit and of the other electrons on a single electron's dynamics. The phenomenological 
constant a describes the relaxation of the self-consistent field [|]. The self-consistent field 
is related to the voltage U across the SSL by E sc = U/l. The current density through 
the SSL consists of two parts: the displacement current jdi sp = ^-E sc (eo is the average 
dielectric constant for the SSL) and the current of ballistic electrons j = —eNV, where V is 
the solution to Eqs (p]). Kirchoff's equation for the resistively shunted SSL, ((e /47r).E sc — 
eNV)S + (E sc l)/R = 0, provides the circuit's contribution to the damping of the self- 
consistent field, a C i r = (RC)^ 1 , C — xf i the total damping rate is the sum of the rates 
from the external circuit and internal mechanisms. Introducing the scaled variables T(t) = 
eaE sc (t)/H, v = rrioVa/h, and w = (8 — A/2)(A/2) _1 , for which the lower (upper) edge 
of the miniband corresponds to w = — 1 (w = +1), we obtain the self-consistent set of 
equations 

v = (J r (t) + ijj s cos Vtt)w — ^ v v, w = — {J-{t) + uj s cos Qt)v — / ~f £ (w — w^), T = uj 2 c v — a J- '. 

(2) 

Here uj s = (eaEo)/Ti and uj c = [(2ne 2 N a 2 A) / \U 2 : eo)] 1/ ' 2 , where uj c is a characteristic frequency 
that can be interpreted as a miniband plasma frequency 0. Throughout this paper we 
assume that initially the electrons are at the bottom of miniband, w(0) = —1, with both 
current and voltage absent, v(0) = ^"(0) = 0, and that the electrons relax to the bottom of 
the miniband, i. e. w {0) = -1, where w (0) = (£ (0) - A/2)(A/2)" 1 . 

Since we are here interested primarily in stationary solutions of @, we do not consider 
transient effects. The appearance of a dc current and a dc bias indicates the presence of 
the zero- th Fourier harmonic in the stationary solutions for v(t) and J-{t). Hence, we first 
establish how the existence of these zero-*' 1 harmonics relates to the symmetry properties of 
(0). Note that the appearance of these zero-* h harmonics is accompanied by other (normally 
forbidden) even harmonics, e.g. the second harmonic. 

The following transformation leaves the system (^|) invariant: 

v-^—v, T -> - J 7 , t^t + T sym /2, T sym = (2n + 1)T, (3) 

where n is an integer and T = 2it/VL is the period of the external ac field. In other words, for 
any solution of system (§) {v(t), w(t), J~{t)}, there exists a symmetry-related "conjugate" 
solution {— v(t + T sym /2),w(t + T sym /2), — T{t + T sym /2)}. Physically, transformation ([3]) 
demonstrates that both signs of electron velocity v and self-consistent field (voltage) T are 
equivalent in an SSL driven by an ac field without bias; in contrast, the lower and upper 
edges of the miniband are not equivalent, thus w -/-> —w. 

Importantly, the system (FJ) possesses two well-studied limiting cases 0. First, when 
7„ = 7 £ = 0, it can be transformed to the resistively shunted junction (RSJ) model familiar 
from the theory of Josephson junctions |]10|| . Second, in the absence of an applied field 
(oj s = 0), but for arbitrary damping, it can be shown || to be equivalent to the familiar 
Lorenz model [11]. In both these limiting cases, symmetry properties have been applied 
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effectively to characterize and interpret solutions, and we therefore expect this to prove true 
in the present, more general case. 

The form of the symmetry transformation suggests the use of a stroboscopic map of 
period T/2 to characterize the solutions of system (||]): This amounts to considering v, w, 
T at discrete times t n = nT/2, where n is an integer. We are especially interested in the 
projection of this stroboscopic phase portrait onto the v — T plane. Let <f>(t) be one of the 
variables v(t) or T(t) and denote by ((f)) the value of <f> averaged over a large time interval. 
By analogy to the RSJ model [|i~0| , |T2|| and the Lorenz model with periodic external driving 
[TT| , we have the following natural classification for attractors of (Q) based on their symmetry 
properties: 

(i) A symmetric limit cycle is a solution that obeys the equality —<f>(t + T sym /2) = <f>(t), 
which implies that it is periodic with period T sym = (2n+l)T; as this period is (2n+l) times 
the external driving period, this limit cycle is also known as a period (2n+ resolution. The 
T/2-stroboscopic plot consists of 2(2n+ 1) points, and its projection on the v — T plane has 
inversion symmetry about the origin. More importantly, ((f)) = 0, and therefore dc current 
and bias are absent. The simplest example of a symmetric limit cycle is a period-1 solution 
of the form <f)(t) = cos Qt + B$ sin Qt, where and are constants. 

(ii) For a symmetry-broken limit cycle, the T/2-stroboscopic phase portrait in the v — T 
plane lacks inversion symmetry about the origin. The simplest example is a period-1 solution 
of the form <f>(t) = A^ cos Qt + B^ sin Qt + C^/2 with zero harmonic equal to (A^, B^ and 
Ctj, are constants). This solution satisfies the equality —(fi(t + T sym /2) = (f>(t) — C^. Actually, 
there are two symmetry-related conjugate limit cycles differing from each other only by 
the sign of the constant C^, and for fixed parameters, a trajectory selects one limit cycle 
depending on its initial conditions. For our case with fixed initial conditions, a trajectory 
can switch to a limit cycle with different symmetry only when the parameters are varied. 
This means that for symmetry-broken limit cycles dc current and bias are nonzero and their 
signs are fixed for fixed values of the system parameters. 

(iii) For symmetric chaos, the solutions are aperiodic, and the T/2-stroboscopic phase 
portrait has a fractal structure with an infinite number of points but possesses inversion 
symmetry about the origin of coordinates in the v — T plane. Such solutions correspond to 
strange attractors with (<f>) = 0, and no dc current. 

(iv) For asymmetric chaos, the fractal pattern in the T/2 stroboscopic representation lacks 
this v — T inversion symmetry, and one has a symmetry-related conjugate pair of strange 
attractors. For asymmetric chaos, (<f>) ^ 0, and dc current is spontaneously generated. 

We observe all four types of attractors in our model. The most unusual type of behavior 
- asymmetric chaos - is illustrated in Fig. |I|, which shows a stroboscopic phase portrait of 
an attractor lacking inversion symmetry about the origin. The symmetry-related conjugate 
attractor exists but is not depicted in this Figure. 

The above symmetry-based classification of attractors demonstrates that stationary so- 
lutions both with and without dc current can be realized in our model. To determine where 
in parameter space the different types of attractors exist, we have undertaken numerical 
simulation of system (Q) holding the dissipation rates fixed but varying the external field 
strength and frequency. For modern SSLs, the velocity relaxation times range from tens 
of femtoseconds up to picoseconds. If we take 7" 1 = 0.35 ps 0, then at A = 22 meV, 
a = 42 nm, N = 1.2 X 10 16 cm -3 and e ~ 13 (GaAs), we have ^ v /uj c = 0.1. As a rule, the 
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relaxation time for energy is an order of magnitude longer than for velocity fl3 |, therefore 
we take 7 e /co> c = 0.01. We choose a/uo c = 0.01, which corresponds to a resistance of R = 5.4 
kOhm and an SSL self-capacitance of C = 0.64 fF (S = 10" 7 cm 2 , I = 40 x a = 1.68 x 10" 4 
cm). For these parameters, the locations of broken-symmetry attractors are shown in the 
Q-u s plane in Fig. |2|. To distinguish genuine symmetry-breaking from long transients and 
rounding errors, we take advantage of the observed "phase- locking" and mark as symmetry- 
broken those solutions having | (J 7 ) \> Q/A. There exists a small region of asymmetric 
chaos with (v) and (J 7 ) ^ 0. However, the majority of the symmetry-broken attractors are 
symmetry-broken limit cycles with higher even and odd harmonics as well as the zero- th 



harmonic [T^]. We also observe some solutions with subharmonics, corresponding to period 
doubling. 

The dependence of the spontaneously generated bias and dc current on the frequency of 
the external field is indicated in Fig. || Recall that (J 7 ) is the long-time average of the scaled 
self-consistent field, i.e. it reflects the dc component of the spontaneously generated electric 
field in the SSL and is in fact equal to the "induced" Bloch frequency (u>g = eaE dc /h = (J 7 )). 
We observe from Fig. |3] that (J 7 ) is approximately equal to an integer multiple of Q. This 
"mode-locking" is the analog of the well-known zero-bias ac Josephson steps in a Josephson 



junction |ID,I3|. But in contrast to the RSJ model, the "mode-locking" in our case is only 
approximate - the exact value of (J 7 ) depends weakly on the field strength u s . We observe 

"mode- locking" for periodic limit cycles as well as for asymmetric chaos. The efficiency of 
ac to dc field conversion, i.e. the ratio ujb/uj s =\ Ed c / E$ |, is often in the range 0.5-0.67, and 
the maximal efficiency we observed is ~ 0.9. Therefore, a resistively shunted semiconductor 
superlattice could be an effective THz ac field to dc current converter. 

From (0) follows that for periodic solutions (v ) = auj~ 2 (T). Therefore, if the bias is 

"phase-locked" (J 7 ) ~ nQ (n is an integer), then the dc current is also proportional to 
nQ. Although this formula does not apply to aperiodic motion, we observed the same 
dependence for asymmetric chaos, at least during our simulation time of 17 x 10 3 uj^ 1 . To 
understand the origin of this "phase-locking" , we recall that the wave functions of a time- 
periodic Hamiltonian with period T can be written as Floquet functions, where ip(x,t) = 
exp(—zet/h)u^(x,t), e is the quasienergy, and u$ is periodic in time with period T flTB] . 



Further, as shown by Krieger and Iafrate [17], when an electric field is applied to a perfect 



lattice, the Bloch wave functions evolve as Houston functions until tunneling to higher 
bands occurs. Therefore, in our single-band model in the absence of scattering the electrons 
would evolve as Houston functions indefinitely. But Zak has demonstrated [TTJ] that for a 
Houston function to be a Floquet state, the dc component of the electric field must satisfy 
the relation ei^ c a = 2imh/T for an integer n. In scaled units, this equation becomes 
precisely the "phase- locking" condition (J 7 ) = nQ; intuitively, this corresponds to resonant 
photon absorption between Stark ladder levels separated by an integer multiple of the photon 
energy. Our results demonstrate that this picture generally continues to be valid with the 
inclusion of scattering, when the Stark levels broaden and the "phase-locking" condition 
is only approximately satisfied. In our situation of an SSL without external bias, the self- 
consistent field leads to the spontaneous formation of a Stark ladder with suitable spacing. 

We conclude with two comments. First, we can estimate the size of the spontaneously 
generated current corresponding to present experimental parameters. For Q = 5.7 THz 
{Q/ui c ~ 0.2) and for the typical SSL parameters and dissipation constants mentioned above, 



5 



we find that I dc = (eNSAa/2h)(v) ~ 28/iA ((v) « 0.002) at field strength £ > £j cr) « 4 
kV/cm (&( cr '/u c ~ 1), which should be within the parameter range of current experiments 
J§-[7|]. Second, we should distinguish our results from a recent study by Ignatov et al 



who pointed out that absolute negative conductance in the volt age- current characteristic of 
an SSL irradiated by an ac field of frequency Q could cause an instability of the zero bias- 
zero current state, resulting in the switching to a new state with dc voltage per SSL period 
close to hfl/e. These studies apply to a biased SSL operating in the high frequency regime, 
in which there can be no chaos. In contrast, our studies apply to an unbiased SSL at lower 
frequencies and lower amplitudes of the ac field and in a regime where chaotic behavior may 
occur. 
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FIGURES 



FIG. 1. Projection of the half-period stroboscopic phase portrait on the velocity - self-consistent 
field plane for asymmetric chaos at VL/u c = 0.47, uj s /uj c = 0.78, 7„/oj c = 0.1, 7 e /^ c = cy/uj c = 0.01. 

FIG. 2. A plot of regions of symmetry-broken motion in the external field frequency(f2)-external 
field strength(cj s ) plane. 

FIG. 3. Dependence of the spontaneously generated bias (T) /u c on the external field frequency 
£l/tu c at uj s /uj c = 1.5, j v /lOc = 0.1, 7 e /w c = a/uj c = 0.01 
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